On the Charge Ordering Observed by Recent STM Experiments 
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We present a perspective on the recent discoveries of possible charge ordering in underdoped 
Bi2Sr2CaCu20s+s and NaxCa2-xCu02Cl2 by high-resolution scanning tunnelling spectroscopy. 



In an STM study of optimally doped 
Bi2Sr2CaCu20s+s {Bi2212) at finite applied mag- 
netic field, Hoffman et al discovered that in the vicinity 
of vortex cores the STM tunnelling conductance exhibits 
a checkerboard spatial dependence[l], with Fourier peaks 
at Q = (±27r/A,0),(0,±27r/A) and A « 4.3a (a is the 
lattice constant of the copper oxide plane) . Subsequently 
Howald et al studied optimally doped Bi2212 below Tc 
at zero field, and reported observation of a nondispersive, 
directionally asymmetric tunneling conductance with 
Fourier peaks at Q = (±27r/A, 0), (0, ±27r/A) and A « 4a 
for low energies {\V\ < 25meV)[2]. However this result 
was not confirmed by later experiments [3-5]. In Ref. [3- 
5] high resolution STM studies of near-optimally doped 
Bi2212 below Tc revealed rich bias-dependent quasiperi- 
odic modulations of the STM tunnelling conductance. 
This phenomenon was explained as the quantum inter- 
ference of quasiparticle de Broglie waves[3, 4, 6]. Ref. [7] 
emphasized the deviation in the dispersion of the Fourier 
peaks in the (1,0) and (0,1) crystallographic directions 
from the simplest quasiparticle interference predictions. 

In Ref. [5] Vershinin et al have also reported results 
for T > Tc. Interestingly, at low energies they observed 
bias-independent conductance modulations with Fourier 
peaks at Q = (±27r/A, 0), (0, ±27r/A) and A w 4.7a[5]. 

Two very recent low temperature STM experiments 
on underdoped Sz2212[8] and NaxCa2-xCu02Cl2 
{NaCCOC)]^] reported bias independent quasiperiodic 
modulations of the local conductance over a wide range 
of energy. These modulations have Fourier peaks at 
Q = (±27r/A,0),(0,±27r/A), with A « 4.5a for 5^2212, 
and A « 4a for NaCCOC. Inspection of the real space 
modulation patterns shows that they are in phase upon 
bias reversal[5, 8, 9]. Ref.[10] and Ref.[ll] suggest that 
this bias symmetry contains vital information about the 
origin of these phenomena. Here we would like to point 
out that this symmetry could appear simply because the 
system is a doped Mott insulator. Indeed, doping a Mott 
insulator quite commonly transfers spectral weight (for 
both electrons and holes) proportional to the doping den- 
sity from above the charge gap to low energy. As a re- 
sult, when the doping density modulates in space, both 
the electron and hole spectral functions modulate with 
the same phase. 

Despite the similarities between the phenomena re- 



FIG. 1: a) Caricature of the observed spatial modulation in 
dl/dV in Ref. [9]. b) Spatially averaged spectral function of 
the hole crystal. The spectral peaks marked by arrows corre- 
spond to states whose wavefunction is peaked at the solitons. 



ported in Ref. [8] and Ref. [9], the Bi22l2 and NaCCOC 
results are different in three important aspects. 1) The 
charge period for _Bi2212 is not an integral multiple of 
the underlying lattice constant, while that of NaCCOC 
is. 2) For Bi2212 the 4.5a modulation appears in dl/dV 
spectrum only when \V\ > 65mV, while in NaCCOC the 
signal appears right from zero bias. 3) For \V\ < 30mV 
Bi2212 exhibits bias-dependent quasiparticle interfer- 
ence modulation[3], while such a signal is difficult to dis- 
cern in NaCCOC. Due to these differences it is possible 
that the phenomenon observed in NaCCOC is due to 
true charge order, while that observed in Bi22\2 is due 
to a yet-uncondensed charge order [8, 9]. 

The above STM experiments constitute an important 
set of observations concerning the electronic structure of 
the cuprates. The purpose of this paper is to present our 
interpretation of the results in Ref.[l, 5, 8, 9]. 

For NaCCOC , the dl/dV modulation patterns con- 
sist of clearly visible 4x4 unit cells over a wide range 
of bias voltages [9] {-150mV < V < l50mV). A cari- 
cature of the real space pattern seen in Ref. [9] is shown 
in Fig. la, where bright/dark marks regions of high/low 
differential conductance. In Ref. [9] it is proposed that 
this spatial structure is due to the formation of a Wigner 
crystal of doped holes. In the following we corroborate 
this proposal with a theoretical analysis. 

We start from a very underdoped situation where 
X = 1/16. This is the ideal doping density for a 4 x 4 
unit cell hole crystal [12]. We interpret the bright region 
in each unit cell as having a high hole density (from a 
single hole dclocalizcd in that region), and the dark re- 
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gion as having a low hole density. In the literature it 
has been proposed that doping a Mott insulator creates 
solitons (or spin bags/composite quasiparticles) [13-15]. 
We believe that the above bright and dark regions reveal 
the internal structure of each soliton. Recently Chen 
et al proposed that at commensurate doping densities a 
Wigner crystal of hole pairs can form[16]. We believe that 
due to Coulomb repulsion such a state is energetically 
unstable. Moreover, even with a strong short-ranged at- 
tractive interaction to overcome the Coulomb interaction, 
we would expect segregation into hole-rich and hole-poor 
regions rather than crystallization of hole pairs. 

To demonstrate the possibility of a soliton crystal we 
perform a variational study of the usual Hubbard model 
with Coulomb and nearest neighbor exchange interac- 
tions added[17], 

H = Huubhard + X] '^ij'^i^J + J X] ' 
{ij) <ij> 

The variational ansatz we use is the most general 
Slater determinant which allows spatially inhomogeneous 
charge density, spin density, bond-current density, and 
superconducting pairing. Technically this is equivalent 
to factorizing the four-fermion terms in Eq. (1) into all 
possible quadratic combinations, and then solving the re- 
sulting quadratic Hamiltonian self-consistently. The fol- 
lowing results are obtained numerically for an 8 x 8 lattice 
with four holes and periodic boundary conditions. 

It is quite encouraging that such a calculation does in- 
deed yield a crystalline arrangement of the doped holes 
for a reasonable choice of parameters {t = — O.SeF, U = 
6eV,J = 0.3eV,Vc = O.SeV, where Vc is the repul- 
sion energy two holes experience at the nearest neigh- 
bor distance) [18]. We have checked that as long as the 
Coulomb interaction exists, the stability of the soliton 
crystal is not affected by moderate modification of the 
above parameters. In all the cases we have studied the su- 
perconducting pairing and orbital current order are both 
absent when the holes crystallize. On the contrary we 
find that a non-zero antiferromagnetic order coexists with 
the hole crystal. 

In Fig. lb we show the spatially averaged spectral hmc- 
tion. The spectral peaks marked by arrows correspond 
to states whose wavefunctions are peaked at the soli- 
ton positions shown by the hole density (Fig. 2a). The 
norm-squared wavefunctions of these states is reflected in 
the spatial variation of the local spectral function within 
±150mey of the Fermi level, as shown in Fig. 2b. Al- 
though the detailed distribution of the spectral weight 
within the unit cell differs for the positive and negative 
bias, the gross structure (i.e. the fact that the overall 
spectral weight peaks near the soliton) is the same. In 
that regard the spatial structure is similar to that ob- 
served in Ref.[9]. Contrary to experiment, in our calcu- 
lation the dark region basically has the LDOS spectrum 




FIG. 2: a)Holc density distribution of the hole crystal. 
b)Spatial modulation of the local spectral function. The im- 
ages shown here are typical for —150meV < E < 150rneV . 
Although the detailed distribution of the spectral weight 
within the unit cell differs for the positive and negative bias, 
the gross structure (i.e. the fact that the overall spectral 
weight peaks near the soliton) is the same. 

of the undoped insulator. We attribute the 150meV 

peak in the dark region seen in Ref. [9] to the extra holes 
running on top of the Wigner crystal background which 
are absent in our calculation. 

We also note that when a soliton is produced on a 
Neel background, its location (on the A or B sublattice 
of the antiferromagnetic order) depends on the spin of 
the removed electron. For example, if removing a spin 
up electron produces a soliton on the A sublattice, then 
removing a spin down electron will produce a soliton on 
the B sublattice. Moreover, our calculation shows that 
after removing an electron of a given spin, an electron of 
the opposite spin is always deeply bound to the soliton. 
This result is consistent with the notion that a soliton is 
a composite particle of a charged void and a spin[14]. 

A weakness of the above variational study is that it un- 
derestimates the qiiantiim fluctuation of spins. To study 
the internal structure of each soliton more carefully, we 
diagonalize the t-J model on a 4 x 4 plaquette with pe- 
riodic boundary conditions. The Hamiltonian we use for 
this study is 

H=Y1 -tPi{clcj, + h.c.)Pj + J(Si • Sj - \nini).{2) 

<ij> 

In the following we choose t = 0.3eV, J = O.leV. Mo- 
tivated by the experimental flndings we restrict the hole 
to a 3 X 3 plaquette in the unit cell. We compute the 
low-energy positive-bias[19] electron spectral function 

Ota 

Here ji is the chemical potential, |\E'o'' > is the one-hole 
ground state with energy Eq^, and |^'^'' > are the eigen- 
statcs of the half-flUed system, with energy at most 
2 J above the half-filled ground state[20]. Fig. 3 shows 
the hole density and A(x, w) at the three inequivalent 
sites of the central 3x3 plaquette. The hole density 
displays a similar spatial structure to that seen in the 
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FIG. 3: Hole distribution and positive bias local spectral func- 
tion. The center of each plaquette is a Cu lattice site. The 
seven boundary sites of the 4x4 unit cell axe replicated using 
periodic boundary conditions. The zero of the horizontal axis 
can be shifted by the chemical potential /n (here /ti = 0). 

Hubbard model calculation. The tunnelling spectrum at 
the brightest site (i.e. site 3) exhibits two peaks. To 
understand these two peaks we note that in a 4 x 4 sys- 
tem the one-hole ground state is a spin doublet. The 
unpaired spin 1/2 is meant to mimic the dangling spin 
in the soliton. After the tunnelling process the spin from 
the added electron can form a singlet or a triplet with the 
original soliton spin. The lower(higher) peak corresponds 
to the lowest energy spin doublet spin singlet (triplet) 
transitions. The peaks seen in Fig. 3 arc similar to those 
observed in the bright region of the 4x4 unit cell in 
Ref.[9]. (Because we constrained the hole to the central 
3x3 region, we cannot address the \^QmeV peak ob- 
served in the dark region of the 4x4 unit cell.) 

In the above discussions we have concentrated on dop- 
ing density x = 1/16, where the system is a charge insu- 
lator. The actual doping density in the sample studied in 
Ref.[9] is likely to be higher [21]. We assert that when x 
deviates from 1/16 the extra holes are dclocalized on top 
of the Wigner crystal background. These mobile holes 
metallicize the system, and at low enough temperatures 
such a metal can become a superconductor. When that 
happens the Wigner crystal order coexists with supercon- 
ductivity. Clearly an important question is the range of 
stability (in doping) of the Wigner crystal. Currently we 
do not have an explicit answer to this question. However, 
we emphasize that in answering this question the effects 
of lattice pinning (including the effect of periodic poten- 
tial as well as elastic lattice relaxation) are extremely 
important. We propose that the 4x4 unit cell Wigner 
crystal with extra mobile holes is adiabatically connected 
to the system with partially gapped Fermi surface near 
the antinodes discussed later in this paper [22]. 

In the following we interpret the STM results for 
Bi2212 [1, 3, 5, 8]. It is widely accepted that at opti- 
mal doping and lowest energies the only important ele- 
mentary excitations are nodal quasiparticles. (The phase 
fluctuations are gapped due to the 3D Coulomb inter- 
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FIG. 4: The non- linear processes by which a quasiparticle is 
scattered by the Friedel oscillation set up by an impurity. 



action.) As the excitation energy increases, antinodal 

quasiparticles and vortcx-antivortex excitations join the 
list. The smallest vortex-antivortex pair is a roton, which 
appears as a (damped) pole of the dynamic density- 
density correlation function. Because the signatures of 
charge ordering at wavevectors Q « (±7r/A, 0), (0, ±7r/A) 
are observed in the underdoped regime, we assert that 
there exist four discrete roton minima at wavevectors 
Qroton ~ (i^r/A, 0), (0, ±7r/A). These roton minima en- 
hance the DC charge susceptibility at Qroton- We be- 
lieve that this enhancement is due to the particle-hole 
scattering across the nearly nested Fermi surface in the 
antinodal region (with nesting wavevector QnesUng, see 

the left panel of Fig.4), so that Qroton ~ Qnesting- 

Due to this enhancement, Friedel-like modulations in 
the electron density at Qroton are induced by disorder. 
These modulations can in turn scatter the quasiparticles 
and give rise to enhanced quasiparticle interference mod- 
ulations at wavevector Qroton- The lowest order such 
process is illustrated in Fig.4a[6, 7, 10, 23]. We believe 
that the bias-independent local conductance modulations 
in Ref.[l, 5] of optimally doped Bz2212 are due to the 
scattering described above. 

In addition to Fig. 4a the process depicted in Fig. 4b 
can also scatter quasiparticles with momentum trans- 
fer Questing - This proccss is resonantly enhanced if the 
quasiparticle energy is greater than the excitation energy 
of the roton[24]. Given the fact that the A « 4.5a mod- 
ulation phenomenon in Ref.[8] appears for |V| > 65mV, 
we think it quite likely it is due to the scattering pro- 
cess shown in Fig. 4b. Finally, as pointed out by Zhou 
et al[25], the process shown in Fig. 4c (without impurity 
involvement) can be very eS'ective in destroying the co- 
herence of antinodal quasiparticle excitations in ARPES 
studies [26]. 

As a function of decreasing x the roton minima deepen. 
When these minima touch zero, roton condensation oc- 
curs. This marks the onset of charge order. In an 
isotropic superfluid such as ^i?e the onset of crystalliza- 
tion coincides with the loss of superfluidity. However, 
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in a system like the NaCCOC where the charge order 
is commensurate with the underlying crystalline lattice, 
interstitial/ vacancy superconductivity can coexist with 
Wigner crystal order due to lattice pinning effects. In 
such a coexistence phase, the antinodal quasiparticles are 
strongly scattered by Qrotoni hence a gap can open up in 
the antinodal region. However since Clroton docs not con- 
nect the nodes the nodal quasiparticles can remain gap- 
less if the potential from charge order is not too strong. 
Under these conditions the nodal quasiparticle peak can 
survive in the ARPES spectra[27]. Upon further decreas- 
ing X superconductivity eventually goes away. The loss 
of superconductivity is triggered by the Bose condensa- 
tion of single vortices. With single vortices condensed 
all doped holes localize and the system becomes an in- 
sulator. After localization of holes the antiferromagnetic 
order is presumably restored. Because a real system has 
disorder, the Wigner crystal order discussed above can 
only have finite-range correlations due to the Imry-Ma 
physics. 

Is the charge ordering discussed here and in Ref.[5, 8, 9] 

responsible for the opening of the pseudogap? As has 
been pointed out in Ref.[8] this crucially depends the re- 
lation between the charge freezing (crossover) temper- 
ature Tcharge and the pseudogap temperature T*. If 
Tcharge = T* & Strong argument can be made that the 
charge order actually drives the pseudogap formation. 
However, in our opinion it is likely that the pseudogap 
is caused by a different mechanism (e.g. spin singlet for- 
mation), and Tcharge < T* . 

Finally, we note that like stripes [28], the hole crystal 
discussed in this paper represents yet another form of 
charge order. We do not have much to say about how the 
system chooses one form or the other. What is important 
is that the cuprates seem to have a tendency to form some 
type of charge order. Whether these charge orders have 
anything to do with pairing is entirely unclear to us. 
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